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ABSTRACT 

Cygnus X-l (Cyg X-l) is the archetypal black hole (BH) binary system in our Galaxy. We report 
the main results of an extensive search for transient gamma-ray emission from Cygnus X-l carried 
out in the energy range 100 MeV - 3 GeV by the AGILE satellite, during the period 2007 July - 2009 
October. The total exposure time is about 300 days, during which the source was in the "hard" X-ray 
spectral state. We divided the observing intervals in 2-^4 week periods, and searched for transient 
and persistent emission. We report an episode of significant transient gamma-ray emission detected 
on 2009, October 16 in a position compatible with Cyg X-l optical position. This episode, occurred 
during a hard spectral state of Cyg X-l, shows that a 1-2 day time variable emission above 100 MeV 
can be produced during hard spectral states, having important theoretical implications for current 
Comptonization models for Cyg X-l and other microquasars. Except for this one short timescale 
episode, no significant gamma-ray emission was detected by AGILE. By integrating all available data 
we obtain a 2er upper limit for the total integrated flux of Fy t u.L. — 3 x 10 -8 ph cm~ 2 s _1 in the energy 
range 100 MeV - 3 GeV. We then clearly establish the existence of a spectral cutoff in the energy 
range 1-100 MeV that applies to the typical hard state outside the flaring period and that confirms 
the historically known spectral cutoff above 1 MeV. 

Subject headings: gamma rays: stars — stars: individual (Cygnus X-l) — X-rays: binaries 



1. INTRODUCTION 

Cyg X-l is a binary system (discovered by Bowyer et 
al. 1965) containing a 09.7 lab supergiant star orbit- 
ing (5.6 days of period) around a compact star with 
a mass function of / = 0.23 ± 0.01 A/© (Gies et al. 
2008) and a mass lower limit in the range 6 13 Mq 
(Ziolkowski 2005). Cyg X-l is then the only known 
high-mass black hole (BH) binary system in our Galaxy 
(e.g., Tanaka & Lewin 1995), and attracted consider- 
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able attention since its initial mass rang e determinations 
([Bolton 19721 : IWebster fc Murdin 19721 ). Being among 
the brightest X-ray binaries in our Galaxy (for a rela- 
tively small distance of 2 kpc and average sub-Eddington 
X-ray luminosity for a 10 solar mass compact object), the 
system has been extensively monitored in the radio, IR, 
UV and X-ray energy bands (see Zdziarski & Gierlihski 
2004 for a review). 

The system spends most of its time in the so called 
"hard state" characterized by a relatively low flux of 
soft X-ray photons (1-10 keV), a clear peak of the pho- 
ton energy spectrum in the hard X-ray band (around 
100 keV), and an energy cutoff around 1 MeV (e.g., 
Gierlihski et al. 1997, McConnell et al. 2002, Del Monte 
et al., 2010). Occasionally, Cyg X-l changes state shift- 
ing its energy power spectrum to a "soft state" charac- 
terized by a large flux in soft X-rays, a lower hard X- 
ray flux, and a tail extending to ene rgies up to 1 MeV 
and beyond (jMcConnell et al. 20021) . Cyg X-l is also 
detected in "intermediate hard states", which usually 
show a less intense hard X-ray emission and a shift of 
the spectral hump towards energies le ss than 100 keV 
(jMalzac et al. 20061 : Iwilms et al. 20061 ). 

Variability in Cyg X-l above 100 keV was observed 
on several different time scales, from months to mil- 
liseconds (e.g. Brocksopp et al. 1999, Ling et al. 
1997, Pottschmidt et al. 2003, Zdziarski & Gierlihski 
2004) and giant outburst episodes have been detected 
in the 15-300 keV by the Interplanetary Network 
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( Golcnct skii et al. 20031) during both spectral states. 

Theoretically, accretion processes onto a BH system 
are extensively studied using Cyg X-l as a typical ex- 
ample. In particular, disk hydrodynamics and radia- 
tive and pair-creation properties of Cyg X-l have been 
modeled with particular emphasis on the X-ray range 
and the highest detectable energies (e.g., Zdziarski 1988; 
Gierlihski et al. 1999; Bednarek & Giovannclli 2007; 
Zdziarski, Malzac & Bednarek 2009). Extensive mod- 
elling of Cyg X-l X-ray spectral states have been car- 
ried out using Comptonization models (e.g., Titarchuk 
1994; Poutanen & Svensson 1996; Coppi 1999) and inter- 
pret the historical data available in the literature with a 
spectral cutoff near 1 MeV. Since the detection of a non- 
thermal power law spectral component extending up to 
~ lMeV energies during the "soft" and "intermediate" 
states, the issue of determining the variability and high- 
est photon energies from Cyg X-l has been of crucial 
theoretical importance. A detection of photon emission 
well above a few MeV from Cyg X-l would provide a 
clear signature of efficient non-thermal acceleration pro- 
cesses occurring in the system, that would need to be 
accounted for in Cyg X-l models and BH accretion disk 
modelling. 

Before the AGILE extensive monitoring of Cyg X-l, 
only temporally sparse information has been available 
in the energy range above a few MeV. The gamma-ray 
instruments on board of CGRO observed the Cygnus re- 
gion several times (typically with 2-^4 weeks long in- 
tegrations) during the period 1991-1997. In particu- 
lar, the EGRET instrument provided an overall up- 
per limit to the flux of 10 x 10 -8 phcm~ 2 s _1 above 
100 MeV. EGRET observations occurred always during 
"hard" spectral state, and did not cover at all the "soft" 
state. 

The only observation of CGRO during a soft state of 
Cyg X-l was carried out in Ju ne, 1996, followi ng an X- 
ray alert provided by RXTE dCui et al. 19971) . OSSE 
and COMPTEL observed Cyg X-l from June 14 to June 
25, 1997 and this led for the first time to the detec- 
tion of a high energy power-law up to about 7 MeV 
(jMcConnell et al. 20 02). This indication of a power-law 
component extending to MeV and beyond was also sup- 
ported in recent years by several INTEG RAL observa- 
tions of Cyg X-l (iCadolle Bel et al. 20061) . 

A remarkable, although isolated, TeV flaring event of 
very high-energy emission above ~300 GeV from Cyg X-l 
was reported by the MA GIG C herenkov telescope during 
a set of observations in 2006 (jAlbert et al. 20071 ). The 
reported VHE emission (for a pre-trial significance above 
4cr) was detected on 2006, September 24, for about 1 
hour (corresponding to an orbital phase of 0.9) during 
a relatively bright hard X-ray emission phas e. Simulta- 
neous INTEGRAL data (jMalzac et al. 20081 ) show that 
the TeV flare from Cyg X-l was detected ~ lday before 
an intense peak in hard X-rays. However, at the time 
of the TeV flare, both the soft and hard X-ray emission 
do not show significant variations or rapid state changes: 
the spectral state was a "hard" one. This detection of 
transient and very rapid TeV emission from Cyg X-l in- 
dicates that extreme particle acceleration processes may 
occur also during a hard spectral state, paving the way to 
detect non-thermal components also in states previosly 
believed to be characterized by a cutoff above a few MeV. 



In this Letter we report the AGILE search for short 
(days- weeks) timescale gamma-ray emission from Cyg X- 
1 in the energy range 100 MeV - 3 GeV with a total 
exposure time of ~ 300 days, during the period 2007 
July - 2009 mid-October. Our data provide the first 
long timescale monitoring for this important BH system. 
A separate paper (Del Monte et al. 2009) addresses the 
details of the X-ray emission as monitored by AGILE 
and other detectors during our first year of observations. 

2. AGILE 2007-2009 OBSERVATIONS OF CYGNUS X-l AND 
DATA ANALYSIS 

The AGILE mission h as been operating since 2007 
April (jTavani et al. 2008| ). The AGILE scientific in- 
strument is very compact and is characterized by two 
co-aligned imaging detectors operating in the energy 
ranges 30 MeV - 30 GeV (GRID, Barbiellini et al. 
2002, Prest et al. 2003) and 18-60 keV (Super- 
AGILE, Feroci et al. 2007), as wel l as by an antico- 
incidence system dPerotti et al. 20061) and a calorimeter 
(jLabanti et al. 20061 ). AGILE 's performance is charac- 
terized by large fields of view (2.5 and 1 sr for the gamma- 
ray and hard X-ray bands, respectively) and optimal an- 
gular resolution (PSF=3° at 100 MeV and PSF=1.5° 
at 400 MeV). Flux sensitivity for a typical 1-week ob- 
serving period can reach the level of several tens of 
10^ 8 phcm" 2 s _1 above 100 MeV, and 10-20 mCrab in 
the 18-60 keV range depending on off-axis angles and 
pointing directions (see Tavani et al. 2008 for details 
about the mission and main instrument performance). 

The A GILE satellite repeatedly pointed at the Cygnus 
region for a total of ~ 315 days (~ 13 Msec net exposure) 
during the period 2007 July - 2009 mid-October. The 
analysis of gamma-ray data presented in this paper was 
carried out with the AGILE-GRID FT3ab2Buildl8 cali- 
brated filter with a gamma-ray event selection that takes 
into account South Atlantic Anomaly event cuts and 80° 
Earth albedo filtering. Throughout the paper, statisti- 
cal significance assessment and source flux determina- 
tion was established using the stand ard AGILE mult i- 
source likelihood analysis software (|Chen et al. 20101) . 
The method provides an assessment of the statistical sig- 
nificance in terms of a Test Statistic (TS) defined as in 
Mattox et al. 1996 and asymptotically distributed as a 
X 2 /2 for 3 degrees of freedom (x 2 /2). 

2.1. Search for persistent gamma-ray emission 

Multi-source likelihood analysis was used to search 
for persistent emission from Cyg X-l position in the 
integrated sky map of the Cygnus region above 100 
MeV for the period 2007 July - 2009 October (Fig □] 
upper panel). The region is characterised by AGILE 
gamma-rays data showing two most prominent sources 
1AGL J2022+4032 and 1AGL J2021+3652 detected with 
high confidence (38.8er and 24.6<r respectively) and a 
gamma-ray flux F 1 = (123 ± 4) x 10~ 8 phcm~ 2 s _1 
and F~, = (57 ± 4) x 10~ 8 phcm~ 2 s^ 1 respectively 
(jPittori et al. 20091 ). We also detect Cygnus X-3 (3.2<r, 
F 1 = (10 ± 3) x 10- 8 phcm" 2 s^ 1 ; Tavani et al, 2009), 
and the nearby pulsar source 1AGL J2032+4102 (10.8a, 
F 1 = (35±3) x 10" 8 phcm~ 2 s" 1 ) 21 . No statistically sig- 

21 A GILE flux va lues are in agreement with the Fermi detections 
HAbdo et al. 2009al) for common sources. 
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nificant gamma ray source is detected at a position con- 
sistent with that of Cyg X-l. The 2-sigma upper limit 
for the gamma-ray flux in the energy range 100 MeV - 3 
GeV is equal to 3 xl0 _8 phc:m~ 2 s~ x . 

Data integrations of 2-^4 weeks exposure from single 
observation blocks give typical 2-a upper limits in the 
range (10-30) x 10" 8 phcm~ 2 s -1 . 

2.2. Search for transient gamma-ray emission 

Motivated by the X-ray variability of Cyg X-l and by 
the particular seq uence of flaring ga mma-ray emission 
from Cygnus X-3 (|Tavani et al. 2 009). we carried out a 
systematic search for short (day) timescale variability of 
the Cyg X-l gamma-ray emission. We used two indepen- 
dent and automatic methods for a blind search of can- 
didate gamma-ray transients in the region surrounding 
Cyg X-l. 

1. The AGILE-GRID multi-source Likelihood method. 
The standard analysis pipeline uses a multiple source 
likelihood analysis that iteratively optimizes position, 
flux and significance of each source by successive rep- 
etitions in which the parameters of one source are varied 
keeping all the others fixed. This method is very effi- 
cient for relatively strong sources and takes into account 
the Galactic diffuse e mission and residual background 
(|Bulgarelli et al. 20081 ). It provides a pre-trial assess- 
ment of statistical significance that needs to be corrected 
when used in repeated systematic searches. For this rea- 
son we also developed an independent method, that takes 
into account multiple comparison corrections (see below). 

2. The False Discovery Rate Method (FDRM). We de- 
veloped a detection method based on the False Discovery 
Rate technique (FDR, Benjamini et al. 1995; Miller et 
al. 2001, Hopkins et al. 2002) that is a statistical test 
taking into account the corrections for multiple testing, 
as needed for example in repeated systematic searches. 
The FDRM allows to control the expected rate of false 
detections (due to background fluctuations) within a se- 
lected sample. The method was adapted to the anal- 
ysis of AGILE gamm a-ray data of the Galactic plane 
(|Sabatini et al. 20101 ) . Given an observed distribution 
of background counts-per-pixel (the null hypothesis), the 
selection is based on choosing pixels characterized by p- 
values 22 smaller than a threshold, cifdr- The crucial 
FDRM feature is that a p-value threshold is not fixed 
a priori (as in traditional statistical methods), but is 
estimated on the data with the requirement that the 
rate of false detections, within the selected sample, is 
the chosen ctF DR. or smaller. A typical value used i n 
the literature (|Miller et al. 2001t iHopkins et al. 20021 ) . 
and that we adopt as a starting value for our search, 
is aprjR=0.05. The FDRM ensures to control this rate, 
while accounting for the ' post-trial' correction of a sin - 
gle detection significance ([Benjamini fc Hochberg 19951) . 
We apply the FDRM in two different ways. 

The global-FDRM (G-FDRM): in this case we carry out a 
blind search for (persistent or transient) sources in large 
(global) daily counts maps of the Galactic plane (0.5° 

22 Given a statistical distribution, a "p- value" assigned to a given 
value of a random variable is denned as the probability, when the 
null hypothesis is true, of obtaining that value or larger. 
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Fig. 1. — AGILE gamma-ray intensity maps above 100 MeV 
of the Cygnus Region in Galactic coordinates displayed with a 3- 
bin Gaussian smoothing. Upper panel: AGILE 2-years integrated 
map. Pixel size is 0.1°. We overlayed the nominal position of Cyg 
X-l (white circle) and the other sources from AGILE catalogue. 
The color bar scale is in units of photons cm -2 s _1 pixel -1 . Lower 
panel: AGILE 1-day map of the flaring episode of Cyg X-l (2009- 
10-15 UTC 23:13:36 to 2009-10-16 UTC 23:02:24). Pixel size is 
0.5°. The black circle is the optical position of Cyg X-l and the 
green contour is the AGILE 2a confidence level. 

pixel size). The null hypothesis for these daily maps is 
the (background dominated) counts distribution of the 
Galactic plane (|6| < 5°). The random fluctuations and 
the diffuse gamma-ray emission of these daily maps are 
well described by Poissonian distributions in AGILE- 
GRID data. Candidate sources in the daily maps are 
identified as significant deviations from the average dis- 
tribution that applies to that specific day. In our analysis 
we use a threshold of oifdr = 0.05 that limits the con- 
tamination by false positive sources in the sample below 
5%. 

The source-FDRM (S-FDRM): the S-FDRM searches for 
flaring episodes in the counts' light curve extracted from 
the position of a single candidate source location. In the 
(verified) assumption that the average source flux at a 
given position is typically below the instrument sensitiv- 
ity, unless it is producing (rare) flares, the null hypothesis 
in this case is obtained by measuring the distribution of 
photon counts for the specific sky location observed at 
intervals of 1-day. We considered the nominal Cyg X- 
1 position and used an aperture search radius of 1.5°. 
As in the case for G-FDRM, candidate flaring sources 
are detected as deviations from the Poissonian average 
distribution, i.e. fluctuations with p— value below the 
chosen threshold. 

3. THE GAMMA-RAY FLARE OF OCTOBER 15-16, 2009 

All of the available AGILE data in the archive from 
2007 June to 2009 mid-October were searched for vari- 
ability on timescales of 1-day with both the Likelihood 
and FDR methods. We used only data within 40° from 
the pointing direction and removed all data affected by 
non-nominal satellite pointings. In this Letter we con- 
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sider only candidates with at least 5a pre-trial signif- 
icance. Only one gamma-ray flaring episode was def- 
initely detected in our thorough search by both inde- 
pendent methods. The bottom panel of Fig [1] shows 
the AGILE gamma-ray intensity map above 100 McV 
for this episode. The emission peaked during the time 
interval 2009 October, 15 (UTC 23:13:36) to 2009 Oc- 
tober, 16 (UTC 23:02:24). The AGILE-GRID multi- 
source likelihood analysis finds a TS=28.09 (=5.3(7 pre- 
trial, 4<7 post-trial 23 , according to xi/2 distribution 
and multiple testing correction) detection at the posi- 
tion (l,b)=71.2, 3.8 ± 0.7 (stat) ± 0.1 (syst) consistent 
with the position of Cyg X-l, for a gamma-ray flux of 
F 1 = (232 ± 66) x 10~ 8 ph cm~ 2 s _1 in the energy range 
100 MeV - 3 GeV. The detection is validated by both 
FDR methods: the G-FDRM analysis finds the source 
with cifdr = 0.05 and the S-FDRM analysis with a 
highly significant cxfdr = 0.001. G-FDRM detection 
has a lower significance due to the use of an average 
background distribution which in this case overestimates 
the local background. For comparison, during the same 
time interval the source 1AGL J2022+4032 (Pittori et 
al. 2009), apparently coincident with the SNR Gamma- 
Cygni, is detected with 3.1<r significance with the likeli- 
hood analysis (and olfdr-g = 0.05), and a flux of F 7 = 
(155 ±60) x 10- 8 phcm- 2 s _1 . The Super- AGILE (18-60 
keV) flux for Cyg X-l for the day is F SA = (580 ± 48) 
mCrab and the ASM flux is F ASM = (268 ± 20) mCrab 
in the 2-12 keV range. The spectral state of the source 
was determined by means of the colour-colour diagram 
obtained from ASM data as discussed in Del Monte et 
al. 2010. Interestingly the flaring episode (MJD= 55120) 
occurred during a hard spectral state. The orbital phase 
of Cyg X-l was in the range 0.38-0.56. The system was 
detected to subsequently evolve into one of the relatively 
rare dips of the hard X-ray light curve. 

4. DISCUSSION 

Fig [5] shows the X-ray historical light curves of Cyg X- 
1 from 2005 November: the upper panel reports the 
RXTE/ASM data, and the lower panel the Swift/BAT 
data, superimposed with the Super-AGILE data from 
2007 November (gray dots). Gray zones highlight AGILE 
pointings of the Cygnus region. The ASM data, show that 
after MJD 53900 the system did not undergo clear tran- 
sitions to one of its soft states anymore. The Swift/BAT 
hard X-ray data are available for the last 4 years and 
show a pattern with rare dips occurring almost once a 
year. 

The AGILE dataset extends for ^300 days, during 
which the system was in its typical hard X-ray state 
(|Del Monte et al. 2010T ). The lack of relatively strong 
gamma-ray emission on a timescale of weeks together 
with the deep upper limit obtained by integrating all 
AGILE-GRID data clearly confirms the existence of a 
spectral cutoff between 1 and 100 MeV in the typical 
hard state. Fig. [3] (upper panel) shows the spectral en- 
ergy distribution of Cyg X-l with its typical historical 
spectral states In the same figure, typical AGILE upper 
limits are given for 2-, 4- weeks and ~ 300 days inte- 
grations. This gamma-ray average spectral behaviour of 




23 This corresponds to a p- 
5.2 X 1CT 4 (post-trial). 



value of 1.7 X 10 6 (pre-trial) and 



Fig. 2.— Upper panel: RXTE/ASM daily light curve of Cyg X- 
1 over the period 2005 November to 2009 October, in the energy 
range 2-10 keV. The black arrows show the day of MAGIC and 
AGILE flaring episode. Gray regions are AGILE pointings of the 
Cygnus Region. Lower panel: Swift BAT long term daily light 
curve in the energy range 15-50 keV and Super-AGILE data (gray 
dots) when available. 



Cyg X-l in the hard state during week- month timescales 
is in overall agreement with Comptonization models of 
black hole candidates (e.g., Titarchuk 1994; Poutanen & 
Svensson 1996; Coppi 1999) and more specifically of Cyg 
X-l (Gierlihski et al. 1997; McConnell et al. 2002). 

However, our detection of October 16, 2009 is the first 
reported 1-day gamma-ray flare in the energy range 100 
MeV - 3 GeV from the system during a hard state. 
This shows that physical processes can occasionally be 
more complex than predicted by current models. The 
lower panel of Fig. [3] shows the AGILE-GRID gamma- 
ray detection during such flare together with the spec- 
tral shapes characterizing the different spectral states 
for reference. Efficient particle acceleration occurs also 
in states characterized by the presence of a hot corona 
that should be in pair-Comptonized equilibrium (e.g., 
Zdziarski 1988; Zdziarski et al. 2009). The gamma-ray 
emission can have leptonic or hadronic origin (e.g., Pe- 
rucho & Bosch- Ramon 2008), depending on the model 
as well as on the assumptions on the acceleration site 
(close or far from the inner disk and/or jet). Lack of si- 
multaneous TeV data prevents a more complete spectral 
analysis of the gamma-ray flaring event. We note that 
the TeV spectrum reported by MAGIC occurred also 
during a hard state and having a photon spectral index 
a = 3.2 ± 0.6 ([Albert et al. 2007T ) is in qualitative agree- 
ment with our AGILE spectral detection 24 , even though 
the broad-band spectrum may be complex and have sev- 
eral independent components. A theoretical analysis of 
our results is well beyond the scope of this paper. 

24 The gamnma-ray detection has a flux about a factor three 
above the model of Zdziarski et al. 2009. 
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Fig. 3. — Cyg X-l spectral energy distribution in typical states ("hard" in solid line, "soft" in dotted line, "intermediate" in dashed 
line; Gierlihski et al. 1999, Zdziarski et al. 2002). The dashed line extrapolated from the hard X-ray state is a purely graphical extension 
of the trend suggested by the historical data. Upper panel: AGILE la upper limits above 100 MeV for integration times of: 2 weeks (A), 
4 weeks (B) and ~ 315 days (C). Lower panel: AGILE data above 100 MeV for the flaring episode. 



5. CONCLUSIONS 

AGILE extensive monitoring of Cyg X-l in the en- 
ergy range 100 MeV - 3 GeV during the period 2007 
July - 2009 October confirmed the existence of a spec- 
tral cutoff between 1-100 MeV during the typical hard 
spectral state of the source. However, even in this state, 
Cyg X-l is capable of producing episodes of extreme par- 
ticle acceleration on 1-day timescales. Our first detection 
of a gamma-ray flare above 100 MeV adds to the even 
shorter detection in the TeV range by MAGIC. These 
data have great relevance for a more detailed theoretical 
modeling of pair equilibrium Comptonized coronae and 
non-thermal particle acceleration that may co-exist for 
short timescales of order of hours-days. 

We note that the gamma-ray flaring activity detected 
by AGILE from Cyg X-l during its decreasing trend of 



hard X-ray emission is qualitatively similar (transition to 
a hard x-ray minimum) to what obs erved in the case of 
the other microqu asar Cygnus X-3 (|Tavani et al. 20091: 
lAbdo et al. 200 9b). Whether this behavior is common 
to microquasars and BH accreting systems is a fascinat- 
ing question that will be addressed by future observa- 
tions. 
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